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ABSTRACT 

Steady state accretion discs larger than ~ 0.01 — 0.1 pc are known to be gravitationally 
unstable for the accretion rates needed to explain super-massive black hole (SMBH) 
activity. We propose that SMBH are fed by a succession of mass deposition events 
with randomly directed angular momenta. Because of incomplete angular momentum 
cancellation a warped accretion disc forms in the inner few parsec. The orientation 
of the disc performs a random walk. Deposition of new material promotes SMBH 
accretion at rates much faster than viscous. Observational implications of this picture 
include: (i) lighter accretion discs that can fuel AGN and quasars and yet avoid star 
formation at R 3> 0.1 pc; (ii) star formation inside the disc is not a function of mass 
accretion rate only. It can take place at high or low accretion rates, e.g., when too few 
clouds arrive in the inner region. An example of this might be the central parsec of our 
Galaxy, (iii) The discs can form Compton-thick obscuring structures of ~ parsec size 
as required in AGN unification models; (iv) faster black hole growth resulting from 
misalignment of the disc and the black hole spin in the early Universe; (v) Isotropic 
deposition of SMBH energy and momentum feedback in the galaxy bulge. This may 
help explain the high efficiency with which it seems to be operating in the Universe, 
(vi) No correlation between SMBH activity and the presence of kiloparsec scale bars or 
gaseous discs in galactic bulges; (vii) Bodily collisions between gaseous components of 
merging galaxies facilitate production of gas streams feeding the centre of the combined 
galaxy. Mergers should thus be catalysts of SMBH growth, (viii) Conversely, galaxies 
experiencing fewer mergers are more likely to form massive nuclear star clusters than 
feed their SMBHs. 
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1 INTRODUCTION 

A well known difficulty in fuelling active galactic nu- 
clei (AGN) is that the typical angular momentum of 
gas in the galactic bulge is very large compared with 
that of the last stable orbit around a black hole (e.g., 
iKrolik, 199dCombes. 200l|jogee. 20oj ). Assuming that the 
bulge is the reservoir ultimately supplying mass to the nu- 
cleus, the result is presumably a disc with the size of a frac- 
tion of the bulge radius, i.e., a fraction of a kiloparsec. The 
material at the inner edge of the disc would then have to 
give up its angular momentum rapidly enough to be able to 
lower itself in the SMBH potential well. This can perhaps 
be done through the action of st ellar and gaseous bars and 
other gravitational torques (e.e.. IShlosman et al.. 19901 ). 

However, there is no clear observational evidence 
for a link between AGN activity and the presence o f 
"grand design" gas discs or stellar bars l|Combes. 2 003). 
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A further theoretical impasse for large scale gaseous 
discs is that even if material does reach the inner 
parsec scales, theory predicts that these discs are too 
cold and too massive, and thus should be unstable 

to gravitational fragmentation and star formation (e.g., ^_ 

iPaczvnski. 197dKolvkhalov fc Sunvaev. 198dShlosman fc Begelman, 1989fco 

For SMBH feeding, star formation is a disease that threatens 

the very existence of accretion discs outside the "self-gravity 

radius", R > 0.01 - 0.1 pc l|Goodman. 20031) . If gaseous 

discs are turned into stellar discs, SMBHs cannot grow by 

gas accretion. 

Faced with the se and other theoret ical difficulties, 
iGoodman (20031 ) and lKing fc Pringle (20071) suggested that 
AGN many may be fed by a direct deposition of low angu- 
lar material into the region R < i? sg . The main difficulty for 
this suggestion is that the specific angular momentum of the 
orbit at R = R ss is only I ~ 100 pc km/s, much smaller than 
the more typical I ~ 10 4 - 10 s pc km/s at R = 10 2 - 10 3 pc. 

On the other hand, the SMBH mass, though large, is 
only a tiny fraction(~ 0.001) of the stellar mass of the bulge 



2 S. Nayakshin and A. King 



l|Haring fc Rix. 20041 ), and must be an even smaller fraction 
of the original gaseous mass from which the bulge is made. 
It thus seems not impossible that cloud-cloud collisions, su- 
pernova shocks, and galaxy mergers channel such a small 
fraction of bulge gas on nearly radial orbits ending up in- 
side R < _R sg . Unfortunately, a quantitative treatment of 
these processes would be very model-dependent and hardly 
constraining. 

In this paper we note that gas feeding the SMBH does 
not actually need to have such a small angular momentum 
after all. We propose that SMBHs can be fed by clouds with 
angular momentum an order of magnitude higher, i.e. I ~ 
10 3 pc km/s or more, if there is significant cancellation of 
angular momentum by shocks in the SMBH's vicinity. This 
would happen if the angular momentum of the clouds is 
random. Below we show that such a picture predicts the 
formation of a warped accretion disc that can be much larger 
than 7? = i? sg , and we discuss observational implications 
for current observations of AGN and for the SMBH-galaxy 
connection. We feel that this "stochastic cloud" mode of 
SMBH fuelling has numerous advantages over the "grand 
design" models. 



2 CAN GRAND DESIGN ACCRETION FEED 
SMBH? 

iGoodman (20031 ) presented the most complete and up-to- 
date consideration of gravitational stability of different ac- 
cretion disc models. The standard accretion flow model is 
gravitationally unstable beyond the self-gravity radius _R sg , 
which is a function of the accretion rate in the disc, M. 
A rough power-l aw fit to this dependence from Figure 1 in 
IGoodman (20031 ) is 

R sg « 0.01 pc M~ 2/7 , (1) 

where M is in units of Mq year -1 . For astrophysically 
interesting accre tion rates, then, R sg ~ 0.01 — 0.1 pc. 
IGoodman (20031 ) also found that none of the more sophis- 
ticated models avoids becoming self-gravitating at parsec- 
scale distances from the SMBH, unless an unspecified and 
very large energy source is applied to the disc to keep it 
hot enough. The required energy at 10 pc, for example, was 
at least that expected if all of the disc m aterial was repro- 
cessed by massive stars. IGoodman (20031 ) also pointed out 
that even if a source for this energy was found, the disc 
may still be thermally unstable as the cooling times were 
much shorter than the lo cal dynamical time. Furthermore, 
ISirko fc Goodman (20031 ) demonstrated that such a high en- 
ergy liberation rate near an AGN is actually in conflict with 
the observed spectra of typ ical AGN. 

iThompson et al. (2005) proposed a starburst disc 
model which app ears to overcome the difficulties noted by 
IGoodman (2003). These authors suggested parameterising 
the radial inflow velocity at a fraction of the disc sound 
speed, arguing that external torques or spiral density waves 
may deliver the required angular momentum transfer. This 
makes these discs lighter at a given accretion rate. Addi- 
tionally, Thompson et al.'s star formation rate in the disc 
is limited by the action of energy and momentum feedback 
from massive stars and supernovae, allowing some fuel to 
trickle down all the way to the nucleus. 



Without numerical analysis it is difficult to say whether 
this model will work for parsec-scale AGN discs. One prob- 
lem is that of time scales. Gravitational collapse of the disc 
is expected to take place on the dynamical time scale, 

E>3/2 

td ^ = m/2M 1/2 ~ 10 3 i$ 2 M 8 - 1/2 year S . (2) 

Or ' Mbh ' 

IThompson et al. (2 005 ) show that the disc cooling time is 
much shorter than the dynamical time in their model. On 
the other hand, the lifetime of massive stars is at least a 
few million years. Therefore, the supernovae contribution to 
the feedback might be activated too slowly in this model, 
i.e., only when the disc would already have collapsed gravi- 
tationally. 

Feedback due to radiation and outflows from massive 
stars may appear much faster, perhaps as quickly as required 
by equation [2] However it is not obvious that the feedback 
will be spread uniformly enough in the disc. For a disk of 
scale height H, radiation and outflows from a massive star 
would affect a surface area of size ~ irH 2 only. Thus, to heat 
up the whole disc, ~ (R/H) 2 massive stars are required. 
If we assume that one m assive star form s from ~ 500 Mq 
of gaseous material for a ISalpeter (19551 ) IMF, this would 
require disc mass of Mdi sc ~ 5 x 10 6 Mq/i" 2 ,, where h-2 = 
100H/R. The corresponding star formation rate would be 
MdiscAdyn ~ 5 x 10 3 M year' 1 hzlMg /2 R~c /2 . This seems 
to be too high for the central region of only a few parsecs. 

In addition, analytical estimates 

( Navaksh in fc Cuadra. 2 005) and numerical simulations 
l|Navakshin et al.. 20071 ) show that the stellar component 
of star forming discs decouples vertically from the gaseous 
disc when the total stellar mass becomes comparable to the 
gas mass. Since gas cools radiatively and stellar motions 
do not, stellar discs become more vertically extended, as 
in the Galactic disc in Solar neighbourhood. Feedback in 
massive star-dominated discs in AGN would thus occur 
mainly outside the gaseous disc, diminishing the efficiency 
of feedback deposition further. 

Summarising, there does not appear to be a completely 
convincing way to avoid the self-gravity catastrophe for 
parsec-scale massive AGN discs. The fundamental problem 
is that these discs are very inefficient as far as angular mo- 
mentum transfer is concerned. They need to be massive 
enough to provide a high enough rate of SMBH feeding, yet 
they cannot be hot enough to escape fragmentation. The 
simplest conclusion one can draw from this is that while 
such discs are probably very important for transferring huge 
amounts of gas into the inner parts of galaxies, they proba- 
bly fail to fuel SMBHs. 



3 STOCHASTIC SMBH FEEDING 
3.1 Large scale region 

Figure 1 illustrates the schematics of the model we would 
like to explore here. We assume that during the epoch of 
bulge and SMBH growth, the gaseous medium feeding both 
of these can be crudely divided into two components (left 
panel of figure 1). The first is the dominant part of the 
gas with non-negligible angular momentum. This compo- 
nent forms a large scale gaseous disc that mainly forms stars 
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Type 1 



Type 2 



Figure 1. Schematic representation of stochastic cloud AGN feeding. Left: large scale (lOOpc to 1 kpc) distribution of gas. Most of the 
gas has a significant angular momentum and is in a circularised large scale disc that contributes little to AGN feeding. The quasi-spherical 
halo of gas clouds produces clouds on small angular momentum orbits. Right: Zoom-in view of the central region (the inner parsecs). 
Clouds collide in the hatched "collision sphere" , and form a warped disc whose orientation changes stochastically. Lines of view clear of 
the clouds and warped disc yield type 1 AGN classification whereas lines of sight intersecting the disc or clouds present AGN as a type 
2 source. Note that this division into Type 1 and 2 is caused by the inner stochastically oriented disc and is unrelated to the large-scale 
structure of the galaxy. 



Blackhole spin axis 



Rin R ou t 
Outer warped disc 



form of compact cloudfl We reason further that there will 
be continuous production of clouds on nearly radial trajecto- 
ries through cloud-cloud collisions, their gravitational scat- 
tering off each other or off star clusters, supernova driven 
shock waves, etc. We now consider the structure of the in- 
ner accretion flow assuming that these nearly radial clouds 
provide the dominant mechanism of SMBH feeding. 



at 



Figure 2. Illustration of the accretion disc structure in Figure 1. 
The outer disc extends from R ln ~ 0.01 pc to Rout- Accretion is 
driven by cancellation of angular momentum in collisions between 
the disc and infalling clouds. The inner viscous disc joins the outer 
one at R = ijj n . A short viscous time keeps the inner disc fiat and 
aligned with the inner fringes of the outer disc, but not necessarily 
aligned with the black hole spin. The outer disc is warped and 
can be highly disturbed, presenting a large solid angle for AGN 
illumination. 



for the reasons explained in 33 The second component is a 
less massive quasi-spherical one. A key assumption we make 
here is that these two components do not mix completely, or 
else we would only have the gaseous disc. This immediately 
implies that the quasi-spherical component must be in the 



3.2 Collision sphere and disc formation 

The right panel of Figure 1 zooms in on the inner parsecs 
of the model. Let -R co ii be the radius of the sphere cen- 
tred on the SMBH that is "optically thick" to a typical gas 
cloud. Clouds collide, convert their bulk kinetic energy into 
heat which is radiated away rapidly, and become bound to 
the SMBH inside this region. The lower limit on R co n is 
given by the geometrical size of the infalling clouds them- 
selves, as the different sides of the cloud then collide with 
each other. Observed interstellar clouds form structures on 
a ran ge of scales, but scales of few parsec are most com- 
mon (jElmegreen fc Falgarone. 19961 ). We thus expect that 

-Rcoll > 1 pc. 

The angular momentum of the clouds at 7? co ii is at most 
the local circular orbital value, Z co ii = ficoii-Rcoii- As clouds 
collide, they circularise and settle into a disc. The outer ra- 
dius of the disc, R ou t, is smaller than R co \\, and is controlled 
by the degree to which angular momentum is cancelled in 
these collisions. 

To estimate it, define N as the average rate of number of 
clouds entering the collision sphere R = R co ii- The angular 



1 Alternatively, it can be a hot hydrostatically supported gas. 
The gas cooling time must nevertheless be short to feed the SMBH 
at an appreciable rate. This cooling gas would thus end up form- 
ing cool clouds again. 
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momentum of the material entering the sphere is then, by 
random walk arguments, Zf ccd ~ Wl/A^^i where iV feed = 



Suppose that the collision sphere is within the black 



hole sphere of influence, i?bh = GM^/a 2 , where a is the 
bulge velocity dispersion. In that case the outer edge of the 
disc is given by 



Rou 



Rcol\ 
Afccd 



(3) 



Obviously, in the limit of N — > oo, we get R ou t — » 0, as 
expected in the case of a zero angular momentum inflow. 



3.3 Forced accretion in the outer accretion disc 

The AGN disc we consider has a very unusual feeding mode 
as the angular momentum of the material supplied to the 
disc fluctuates with arrival of new clouds. Unlike binary 
Roche lobe overflow systems, the geometry is 3D, with new 
material arriving not necessarily at the outer edge. In partic- 
ular, parts of clouds with particularly low angular momen- 
tum might strike the disc further in than Rout- In a clear 
difference from the standard accretion flow, no angular mo- 
mentum transport is needed in such randomly fed discs to 
promote accretion. Cancellation of the randomly directed 
angular momentum of the new material and that of the disc 
sets up a radial inflow. We call this "forced accretion" as it 
is mediated by external deposition of material. 

Let us estimate the steady-state mass of such a disc. 
The specific angular momentum of a disc containing N equal 
mass clouds that arrived with a random direction of the an- 
gular momentum will be Zdisc ~ Icoii/N 1 ^ 2 , again by random 
walk arguments. If this angular momentum is less than l Bg , 
the disc will move in radially (in this very simple picture) 
to R < R Bg . Thus, N — (/coii/'sg) 2 random cloud deposi- 
tion events are needed to lower the gas to R < R sg region. 
Accordingly, the "stochastic average" mass of the disc will 
be 



Mdisc = MiN ~ Mi 



'coll 



M 1 R C . 

Rse 



(4) 



where Mi is the mass of a cloud. For example, if Mi = 
100 M Q and R coU /Rsg = 300, then Af disc = 3 x 10 4 M Q . The 
accretion time scale is given by the expression 



_ N _ 
N 



-Rcoll 1 

Rsg N 



(•») 



In a quasi steady state, the surface density profile, 
E(_R), is a strongly peaked function for these discs. We have, 
const= M ~ Y}(R)ttR 2 /i acc . As i acc is independent of ra- 
dius, we obtain E(i?) oc R~ 2 . This is a much steeper depen- 
dence than that for discs in which surface density is regu- 
lated by angular momentum transfer, since viscous time typ- 
ically is a strongly increasing function of radius. For exam- 
ple, for a standard gas-dom inated accretion disc, E oc R~ 3 ^ 
ijShakura fc Sunvaev. 19731 ). 

Consider gravitational stability of such discs. In our 
model the heating per unit surface area is Q+ ~ zn 2 R 2 t-^. 
The radiation flux emerging from the disc is <tT 4 [t-|-1/t] _1 , 
where T is disc midplane temperature and r is the optical 
depth of the disc. Since [t + 1/t] > 1 for any t, we can 
estimate the disc temperature to be 



T > 



T,Q 2 R 2 



21 1/4 



Crf a 



(6) 



Further, if we take into account gas pressure only, we ob- 
tain the minimum scale height of the disc as H — c s Q,~ L = 
(fcT/V) 1/2 fr\ Requiring M disk « Ei? 2 > M hh (H/R), we 
arrive at the minimum disc mass required for gravitational 
instability in our model: 



M b 



k 4 R 



1/7 



0.002 



Rn 



M 2 U 



1/7 



(7) 



where M 8 = Mbh/10 Mq, t 5 = t acc /10 5 years. In the lat- 
ter we set fi — m p . Equation [7] allows us to estimate the 
minimum accretion rate at which these discs would become 
self-gravitatin£|: 

(8) 



Meg 
£acc 



Mr. 



year 



Note that the dependence on the radius is quite weak here. 
Apparently, such externally fed accretion flows are able to 
deliver accretion rates of the order of the Eddington accre- 
tion rate of the SMBH if the time scale on which the angular 
momentum of the incoming gas fluctuates is shorter than 



teto ~ 10° years 



(9) 



This time scale can be much shorter than the disc viscous 
time, t yiac ~ a~ 1 (R/ 'Jf) 2 Q~ 1 at p arsec distances from the 
SMBH IShakura fc Sunvaev (1973 ) . 

We can estimate the maximum outer radius of the ex- 
ternally forced discs. Requiring that the time scale for a 
significant angular momentum change i s to to be larger than 
the local dynamical time given by equation [2] we see that 
the disc cannot be larger than about 10 pc. Another way to 
put this result is to say that any disc larger than this will 
necessarily be star forming. 

3.4 Inner viscous disc 

At small radii, the viscous time is shorter than t acc . Within 
that region, the inflow is viscous rather than externally 
driven. The inner radius of the forced discs can thus be es- 
timated as 



n n ni 2 /3 ,4/3 ,2/3 ,,1/3 

Rin ~ 0.01 pc a ' A h_ 2 t 5 ' M 8 ' , 



(10) 



where hi = 100H/R is the aspect ratio of the inner disc. 
This is an estimate only, as hydrogen ionisation instabil- 
ity might in principle be im portant for this region (e.g., 
ISiemiginowska fc Elvis. 1997h and modulate the accretion 
rate onto the SMBH. 

Lodato & Pringle (2007) show that, in a diffusive 
regime, warp propagation occurs on a time scale shorter 
than the viscous time by a factor 02 /a, where «2 is the 
"warp diffusion viscosity coefficient" . They find that 012 has 
a maximum value of a few which is attained for strongly 
warped accretion discs. Hence, in the context of our model, 
the inner region of radial size 

i? flat -0.1pc (^) 2 V 2 3 t 2/3 M 8 1/3 , (11) 



2 a fin ate disc opacity and the radiation contribution to the disc 
pressure, both neglected above, will make this only limit higher 
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is able to flatten out due to local viscous forces. The inner 
flow then consists of the innermost zone (R < Ri u ) where 
the flow is flat and viscosity mediates the inflow, and the 
outer zone (Ri n < R < Ra&t), where the disc is also flat 
but the inflow is driven by the same mechanism as in the 
larger forced disc. The orientation of the inner disc R < Rf\ at 
coincides with the that of the innermost part of the larger 
forced disc. 

Note that Ra^t sets a velocity scale - the Kepler velocity 
at that radius - of v ~ few xlO 3 km s _1 . Since the inner 
disc is flat, it presents a very small solid angle to the AGN 
illumination. On the other hand, the outer disc might be 
much better exposed to the AGN. It is then possible that 
the transition between the inner flat and the outer warped 
disc will be a site for broad AGN optical and UV lines. We 
would then predict, in the context of our model, that the 
broad line features should have velocity widths of order of a 
few thousand km s _1 . 



4 IMPLICATIONS 

Avoiding the self-gravity catastrophe. The picture of 
accretion proposed above suggests that stochastically fed 
AGN accretion discs can extend beyond the self-gravity ra- 
dius _R sg . These discs are not subject to the self-gravity 
catastrophe as long as the gas feeding these has a ran- 
dom sense of angular momentum fluctuating on time scales 
£ < t sto - This then ensures that super- massive black holes 
can grow via gas accre tion, as required by the observations 
i|Yu fc Tremaine. 20021 ). 

Note that ou r model is relat ed to the proposals ma de by 
iGoodman (20031 ) and especially iKing fc Pringle (20071 ) that 
material feeding AGN comes in shots directly impacting the 
self-gravity radius R sg . The latter is very small, 7? sg ~ 0.02 
pc. We therefore think that the picture proposed here is 
more likely. The incoming clouds do not have to carry such 
a small angular momentum, and they first form a disc that 
can be much larger than i? ag . 

Nuclear star formation. In our model, star forma- 
tion in the direct vicinity of a SMBH is a function of 
not only the gaseous mass deposited there but also the 
manner in which that mass arrived there. In particu- 
lar, contrary to steady-state viscous accretion flows, it 
is feasible to have mainly accretion at a high mass de- 
position rate, but mainly star formation at lower rates. 
An example of the latter situation could be the cen- 
tral parsec of our Galaxy, where two young stellar rings 
have appa rently formed in situ about ~ 6 Millio n 
years ago l|Levin fc Beloborodov. 2003]Genzel et al.. 2 003). 
A reasonable explanation of this is an infall of two 
clouds with mass of a few thousand to ~ 10 4 Mq 
l|Navakshin fc Cuadra. 2005lPaumard et al.. 200rj ). Coeval 
infall of several more clouds of this type from random di- 
rections could have led to more angular momentum cancel- 
lation and mainly accretion of gas on Sgr A* instead of star 
formation. 

There is no one-to-one relation between AGN activity 
and nuclear starbursts in this picture, as it is possible in 
principle to have either one separately. On the other hand, 
AGN feeding in our model depends on availability of gas 



clouds on low angular momentum orbits. Star formation 
feedback (outflows and supernovae) may be encouraging 
such orbits via adding large random velocity kicks to gas 
clouds in vicinity. Thus, statistically it is more likely that 
starbursts and AGN would be linked to one another. 

Masing and/or star forming rings in nearby AGN. 

If SMBHs are fed by grand design discs, the surface den- 
sity distribution is normally a continuous well behaved func- 
tion, such as a power-law. In contrast to that, if accretion is 
stochastic, then the disc surface density distribution does 
not have to be a smooth function, particularly at lower 
cloud deposition rates. In particular, rings with radial ex- 
tent A_R ^ R might result. Observationally this might be 
relevant to masing discs in nearby AGN, where the radial 
extent of the emitting region is usually narrow. If rings are 
massive enough, stellar discs with well defined inner and 
outer edges may form in this way. 

Random SMBH jet orient ation, faster early SM BH 
growth. As pointed out by IKing fc Pringle "(2 007). the 
spin of a black hole fed by deposition of clouds with 
randomly directed angular momentum will be frequently 
misaligned with that of the inner accretion disc. This 
leads to a smaller radiative efficiency and a faster black 
hole growth, helping to explain the heavy-weight champion 
SMB Hs observed already a t high redshifts (tbd). Further- 
more, ISchmitt et al. (20021 ) pointed out that jets in radio 
galaxies seem to be oriented (almost) randomly with respect 
to dust discs of these galaxies. This would be natural in our 
model. 

The Mbh-Mb u ig e correlation. Models explaining the cor- 
relation between the observed SMBH masses and the bulge 
masses hosting them as aris ing due to S MBH accretion feed- 
back seem to be promising (|King. 20031 ). These models pos- 
tulate that gas located throughout the bulge can be over- 
heated or swept away by the SMBH feedback. If SMBH were 
fed from a large scale massive gaseous disc with a small H/R 
ratio, it would be very hard to affect that gas reservoir. 

The difficulty in expelling a flat disc is two fold: (i) 
Firstly, the column depth of the gaseous disc with mass M<j 
through the midplane of the disc is ~ Md/{RH). This is 
R/H times higher than the column depth of spherically dis- 
tributed material with same mass, ~ M/R 2 . Hence it is 
harder to affect the disc with feedback of any type, (ii) Sec- 
ondly, the SMBH feedback may well be collimated and ori- 
ented perpendicular to the inner accretion disc. If the inner 
disc is oriented same way as the much larger disc in the 
bulge, then the feedback may miss the disc altogether. A 
numerical illustration of these principles is provided by the 
recent calculation of t he radiation field during build up of a 
high mass protostar bv lKrumholz et al. (20051 ). These simu- 
lations show that if outflow is collimated, then the radiation 
field is also collimated and therefore affects the infalling ma- 
terial significantly less than previously thought. 

Thus, there does not seems to be a way to expel either 
the disc from the bulge or shut off SMBH feeding if this orig- 
inates in a massive flat disc. On the other hand, if material 
feeding the SMBH is distributed quasi-isotropically through- 
out the bulge, and if the direction of the SMBH spin (and 
feedback) randomly fluctuates, then the arguments made by 
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iKing (20031 ) apply and may explain the Mbh-A^buigo corre- 
lation. 

Absence of a correlation between AGN activity and 
presence of bars. If SMBH is not fed by the large scale 
grand design gaseous features such as spirals or bars, then 
the latter have no bearing on AGN activity, as observed. 

Importance of mergers in AGN feeding. Here we have 
argued that SMBH fuelling is driven by infall of gas on 
nearly radial trajectories. For the gas at R ~ i?bui go to as- 
sume such an orbit, the angular momentum of (only some!) 
clouds needs to cancel out almost completely by cloud-cloud 
collisions, for example. A way of putting much more gas on 
such orbits would be to have a major merger, in which the 
gaseous discs of the galaxies collide bodily. Such collisions 
probably generate SMBH feeding rates well in excess of the 
Eddington rate since the gaseous discs can be orders of mag- 
nitude more massive than the SMBHs. 

Nuclear star clusters and SMBH in dwarf galaxies. 

Galaxies going through fewer mergers channel less gas on 
nearly radial orbits. Central black holes in these galaxies 
are thus relatively more "fuel starved" than their cousins 
in galaxies experiencing more mergers. SMBHs in such "no 
mergers" galaxies could then be underweight compared to 
their expected Mbh-^buige or Mbh-c mass. If smaller galax- 
ies go through fewer mergers, as current cosmological simu- 
lations imply, then it is these galaxies, i.e., dwarf spheroidals, 
that are most strongly affected by this argument. The fuel 
stalled in the central region of the galaxy because of in- 
sufficient angular momentum cancellation can be used up 
in nuclear star formation. These ideas might be relevant to 
the recently claimed dichoto my of nuclear star cluste rs and 
SMBHs in observations (e.g.. IWehner fc Harris. 20 06). with 
the former objects present mainly in low mass galaxies, and 
with dwarf spheroidals possibly lacking SMBHs. 

The final parsec problem for SMBH mergers. It is 

found that central black holes in a major merger always 
find their way into the centr al part of the resultin g galaxy 
by dynamical friction (e.g., Beg elman et al.. 19801 ). These 
black holes then form a binary. The SMBH binary contin- 
ues to shrink by expelling stars. However, when the binary 
separation approaches a few parsec, interactions with the 
stellar background become too inefficient. If gas continues 
to pile up in the inner few parsecs of the galaxy through 
stochastic cloud deposition, the SMBH pair can do work on 
this supposedly unlimited gas supply instead of stars, and 
continue to shrink u ntil gravitational r adiation takes over 
(for related ideas see lEscala et al.. 20051 ). 

Unification and obscuration schemes of AGN. 

Warped accretion discs have been claimed to be important 
in the obscuratio n of the inner regions of type II AGNs 
l|Navakshin. 20051 ). Stochastically fed accretion discs are 
generically strongly warped because different mass "shots" 
will likely have not only random angular momentum orien- 
tation and also different circularisation radii. Further, it can 
be shown that time scales for flattening the warps can be 
much longer than t sto (except for the inner disc, see Sj3]4}. 
We have shown that our disc model is gravitationally 



stable for disc masses less than M sg ~ 0.002Mbh (see equa- 
tion [7J. The column depth of a such a disc is then 

E ~ ^ ~ 20 m 8 5/7 R; c 13/7 t; 1/7 , (12) 

Thus, if these discs are strongly warped, they could form 
Compton-thick absorbers with size of at most a few par- 
sec. This agrees with estimates for the obscuring medium 
inferred in AGN. 

Absence of type II LLAGN. In our picture, Low Lumi- 
nosity AGN would deposit mass at a low rate or always with 
the same angular momentum. In either case this would im- 
ply a disc passively circling the SMBH at the circularisation 
radius or perhaps forming stars if the disc becomes too mas- 
sive. In the first case, the disc may well have enough time to 
flatten out by viscous or gravitational torques, whereas in 
the other the disc might be consumed by star formation too 
rapidly to provide a large enough obscuring column depth, 
therefore, in general we predict that LLAGN must be much 
less obscured than brighter AGN such as Seyfert Galaxies. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have considered SMBH feeding in galax- 
ies. We suggest that due the extremely long time scales for 
angular momentum transfer and the loss of gas to star for- 
mation mean that large scale gas discs do not contribute 
directly to SMBH growth. Instead, we argue that SMBHs 
are fuelled by low angular momentu m gas. Our mod el is 
closely related to sugg estions made bv lGoodman (2003! ) and 
IKing fc Pringle f 20071 ). as the core idea is that the angular 
momentum of the incoming gas is small in a time-averaged 
sense. However, our model has additional observational im- 
plications. For example, warped "externally forced" accre- 
tion discs can extend to scales of up to 10 pc (see the end 
of t|3.3p . which is much larger than the self-gravity radius 
R se ~ 0.01 - 0.1 pc. 

One question that we do not address here is whether 
real galaxies feed enough gas on nearly radial trajectories to 
support the growth of their SMBHs. This requires numeri- 
cal simulations of a significant dynamic range and physical 
complexity. 
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